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1. Introduction 

Adopting Weyl's paradigm of a conformally invariant physics [[I]], though without the 
additional Weyl 1-form, Bach put forth an alternative set of fourth-order gravitational 
field equations — the Bach equations 0]. They are derived by varying, on a four- 
dimensional Riemannian space-time (M,g), the unique quadratic conformally invariant 
action that is obtained upon gauging the global conformal group |§, viz. 

s = -\ ( d 4 x^gC abcd c abcd , (l) 



4jm 

where C a bcd is the Weyl tensor.[j] Bach equations are found in differential geometry and in 
mathematical studies of Einstein's equations of general relativity that utilise conformal 
techniques. For instance, Iriondo et al. have recently shown that the integrability 
conditions arising in the null surface formulation of general relativity impose a field 
equation on the local null surfaces which is equivalent to the fulfilment of the Bach 
equations f|. Indeed, the Bach equations constitute a necessary condition for a space 
to be conformal to an Einstein space. 

Weyl conformal gravity — the theory based on the action ([]]) — has attracted much 
interest as a promising candidate for quantum gravity. Beside the attractive feature 
of being the local gauge theory of the conformal group, quantum conformal gravity 
possesses the virtues of renormalisability, due to the lack of scale in the theory, 
and asymptotic freedom. Paradoxically, conformal gravity has not received general 
acceptance since it involves higher derivatives; some solutions of the classical theory 
are expected to have no lower energy bound and, therefore, exhibit instabilities, i.e. 
'runaway solutions' [§]. In other words, ghosts (negative-energy modes) propagate in the 
linearised theory, thereby breaking unitarity. However, beside the fact that the unitary 
issue is still open in the context of theories with higher derivatives, very interesting 
results inferred from the zero-energy theorem indicate that nonperturbative effects 
may alter the spectrum of the asymptotic states, effectively confining the ghosts || ||. 

On the other hand, conformal gravity bears no resemblance to general relativity 
classically: Einstein gravity is not scale invariant and requires the dimensionful Newton 
coupling constant G^. Moreover, being fourth-order, the Bach equations have a larger 
class of solutions than the second-order Einstein equations; those solutions include 
the aforementioned undesirable runaway modes. Notwithstanding this fact, it is a 
remarkable result that Birkhoff's theorem holds in conformal gravity, unlike generic 
quadratic theories |T0|| . Moreover, every Einstein space or every space conformal to an 
Einstein space fulfils the vacuum Bach equations automatically. This property has 
led many authors to claim that conformal gravity and Einstein's general relativity 
are indistinguishable from each other as far as classical tests are concerned: in both 
theories the Schwarzschild-de Sitter line-element is a vacuum solution. However, as 



firstly pointed out by Pechlaner and Sexl [11 and independently by Havas the 



f We adopt Wald's conventions of sign, definitions of curvature tensors (but the extrinsic curvature), 
and abstract index notation unless otherwise stated 101. 
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Schwarzschild-de Sitter metric cannot be the vacuum solution of the Bach equations 
with matter appropriate to the exterior of a positive-mass body like the Sun. (Rather 
than the observed Newtonian 1/r potential, a point-mass M produces a linearly-rising 
gravitational potential.) This casts doubts on the viability of Mannheim and Kazanas's 
conformal gravity |13| (as exemplified by Perlick and Xu's critics (14]]), which is conceived 
as a classical alternative to Einstein's theory and the virtue of which is its ability to 
account for the shapes of galactic rotation curves without invoking dark matter. 

Nevertheless, one may consider that conformal gravity provides one with a better 
approximation than general relativity to a yet unknown quantum theory of gravity, 
before the expected spontaneous breaking of conformal invariance. In that respect, 
it could make sense to examine whether known properties in Einstein's theory also 
pertain to conformal gravity. With regard to cosmology this leads us to consider the 
spatially homogeneous anisotropic spaces of Bianchi type I. They are indeed the simplest 
nontrivial vacuum space-times since Friedmann-Lemaitre-Robertson-Walker (FLRW) 
isotropic cosmological models are conformally flat. In this paper we lay emphasis on the 
mathematical analysis rather than the physical discussion of the solutions: Our aim is to 
determine whether the (Bianchi type I) Bach equations are integrable or not — in a sense 
that will be made precise later — and to find out exact solutions. As far as this system 
is concerned, it turns out that a nice interplay between involutive methods (within an 
appropriate canonical formalism) and techniques that probe the analytic structure of 
the system enables us to go through the integration process and produce the required 
solutions. 

The plan of the paper is as follows. In Section |2| we present briefly the Bach 
gravitational field equations and the associated Bach tensor. We also set our notations 
for the Bianchi type I metric. We devote Section |3| to the canonical formulation and 
involutive treatment of this system. Two of us have already considered the Hamiltonian 
formulation of Bianchi cosmologies in the pure R 2 variant of the general quadratic 
theory ||15|| . Since in this paper we consider the conformally invariant case, we build up 
a Hamiltonian form of the conformally invariant action (|l|) that we further particularise 
to the Bianchi type I system. We then take the ensuing canonical system as our starting 
point for seeking out exact solutions. For that purpose we perform the Dirac Bergmann 
consistency algorithm on specific constraints that lead to an involutive system. In 
Section || taking up the Painleve approach we prove the non-integrability of the Bianchi 
type I system and extract all possible particular solutions, thereby showing that the 
involution algorithm has produced the complete set of exact solutions that can be written 
in closed analytical form. In Section § we characterise our solutions in regard to their 
conformal relationship with Einstein spaces. This paper is a more thorough presentation 
of results reported elsewhere |T6| . 
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2. Conformal gravity and Bach equations 

Variation of the conformally invariant action (|1|) with respect to the metric yields the 
conformally invariant fourth-order Bach equations 

B a b := 2V m V n C ma 6 n + C ma bnR mn = 0. (2) 

In equation @ B a b is called the Bach tensor; it is symmetric, trace-free, and conformally 
invariant of weight —1. Although expression @ has a very compact form this is 
not appropriate for calculational purpose with computer algebra: Computation of 
the contracted double covariant derivative of the Weyl tensor is quite a heavy task 
even for simple metrics. Recently, Tsantilis et al. have provided an algorithm for the 
MathTensor package, which is based on the decomposition of the Riemann tensor in 
its irreducible pieces, and that gives the Bach equations in a much more tractable form, 
especially with regard to cosmological applications flTTj. In accordance with their results 
we write down the following equivalent expression of B a b, 

B a b = — 0(^R a b — —g a bj + -V a Vfe-R + (C m abn + Rmabn + RmbQanj R™™ , (3) 

where the box □ is the d'Alembertian second-order differential operator. Note that the 
equivalence of expressions (|J) and (|3|) can be proved owing to the useful formula 

V m C ma fe c = V ' c L a b — Vf,L ac , (4) 

where the symmetric tensor L ao is defined by 

— (Rg a b - QRab^j , (5) 

and making use of the contracted Bianchi identities. 

The simplest space-times exhibiting nontrivial physical degrees of freedom in the 
conformally invariant gravitational theory based on the action ([!]) are the spatially 
homogeneous anisotropic spaces of Bianchi type I — the isotropic flrw spaces are 
conformally flat. The corresponding Bach equations (||) can be derived with the help 
of symbolic computational packages such as the ExCALC package in Reduce |18|]. The 
first step is to write the metric in such a way that the conformal invariance be manifest 
ah initio. We assume that the Bianchi type I metric be diagonal. We adopt Misner's 
parameterisation and choose the conformal time dt := e _/ Mr rather than the proper 
time r; hence the line element is 



L 



ab 



ds 2 = e 2 " 



where /z, /3 + , and /?_ are functions of time only. The independent Bach equations 
corresponding respectively to the temporal and spatial diagonal components of the Bach 
tensor @ are given explicitly in |Appendix A . 
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3. Hamiltonian treatment of Bianchi type I spaces 

3.1. Hamiltonian formalism and canonical equations 

Unlike what happens in general relativity, where unwanted second-order derivative terms 
of the metric can be discarded from the gravitational action through pure divergences, 
there is no chance whatsoever to weed out such terms in the context of higher-order 
gravity theories. It is fortunate however that a consistent method of building up a 
Hamiltonian formulation of those theories does actually exist; it is a generalisation of 
the classical Ostrogradsky formalism (see, e.g., ||19|| , and references therein). Basically, 
it consists in introducing auxiliary degrees of freedom that encompass each of the 
successive derivative terms higher than first order. Without resorting explicitly to 
this method Boulware worked out a Hamiltonian formulation of quadratic gravity ||20|| . 
Contradistinctively, Buchbinder and Lyakhovich developed a canonical formalism for 
the most general quadratic gravitational Lagrangian in four dimensions, by employing 
the aforementioned generalised Ostrogradsky method [2~I|| . Here we make use of a 
slightly different generalised Ostrograsky construction — as compared to Buchbinder and 
Lyakhovich's formalism — in order to derive a canonical form of the conformally invariant 
action ([1]), which was first obtained by Boulware [2"0|j . 

Assume that space-time be foliated into a family of Cauchy hypersurfaces E t of 
unit normal n a (see, e.g., [|]). The induced metric h a b onto these hypersurfaces is 
defined by the formula h a b = g a b + n a rib. The way the hypersurfaces are embedded into 
space-time is provided by the extrinsic curvature tensor K a b := —h b c V ' c n a = — \C n h a b, 
where C n denotes the Lie derivative along the normal n a (we adopt C n as a generalised 
notion of time differentiation). In a coordinate basis the normal n has the components! 
n a = (1, —N l )/N and n a = (-N, 0) respectively and the metric g a b can be written as 

ds 2 = h ij (dx i + N*dt) ® (dx j + N j dt) - N 2 dt ® dt. (7) 

The quantities N and N l are the standard ADM variables, namely the lapse function 
and shift vector respectively. The 3 + 1 -splitting of space-time enables one to express 
the Lagrangian density £ = — \\f^g C a bcdC abcd only in terms of the quantities defined 



4 

onto Hf For this purpose we make use of the Gauss, Codazzi, and York equations, 

(3) R a bcd = h\h b f h c °h;R k fge + 2K a [d K c]b , (8a) 
h e d hfh g b n a R\ cd = 2V [f K e]g , (86) 
n b n d h e a h f c R abcd = C n K e} + D (e a /} + K eg K 9 f + a e a f , (8 c) 

respectively, where D denotes the covariant derivative onto S t and a is the four- 
acceleration of an observer moving along the normal to T, t , viz. ab := n c V ' c n b- After 
some algebra we obtain the equations: 

1. 
-/ 
2 

% Latin indices i, j, k, . . . refer to spacelike components in S t whereas Greek indices a, f3, 7, . . . refer 
to components with respect to a specific space-time basis, such as the ADM basis with vectors 
e„ := n = (d t — N l di)/N and = 9j. 



h e d hj c h g b n a C a bcd 



h d h f c h g b - -h bd (h eg h f c - h fg h e c ) V [c K d]b , (9a) 
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(96) 



and (|96|) with respect to the ADM basis: 



Krs\t ~ Krt\s j j 



(10 a) 
(106) 



C\ k ^[6 r i S s j 8{--h rt (h ik S s j -h ij 5 s k 

C mn3 = \ [b\b\ - l -h l3 h kl ) (c n K kl + ^ + KK kl + ®R kl } , 

where a stroke ' | ' refers to the components of the covariant derivative in the ADM basis 
and the symbol n indicates a component along e n . Owing to the fact that the Weyl 
tensor identically vanishes in three dimensions, the Lagrangian density in ([!]) reduces to 



£ = -NVh ( 2C mnj C ninj + C nilk C"v 



(11) 



where C ni j k and C nin j are given by equations ( |10dj and ( |10i| ) respectively. 

Now consider that the induced metric hij and the extrinsic curvature Ky be 
independent variables: are introduced as auxiliary Ostrogradsky variables. In order 
to recover the definition Ky = —^C n hij one we must replace the original Lagrangian 
density by a constrained Lagrangian density, 



£ = N^Z + \ ij (C n hij + 2K, 



13)1 



(12) 



with Lagrange multipliers X tJ as additional variables. One must therefore resort to 
Dirac's formalism for constrained systems (see, e.g., ||22|| , and references therein). The 
conjugate momenta are defined as usual, viz. 

dZ 



p J 
Q ij 



dZ 
dZ 



A* J , (13a) 
---2VhC ninj , (136) 

0. (13c) 

(Observe that the variable Q %3 is related to the electric part of the Weyl tensor.) Hence 
there are primary constraints: 

0« = p « - A iJ « 0, $ = InjQ 1 ~ 0, « 0. (14) 

(The second constraint arises because only the traceless part of the 'velocities', (C n Kij) T , 
can be extracted from equation (|10fl) .) Performing a Legendre transformation on the 
extended Lagrangian density fll2"|) we obtain the canonical Hamiltonian density 

53 C = -2^ + VhC mjk C^ k - -^=Q^Ql - Q T ^.. - Q T ^% - KKjjQ Tij , (15) 

which is vanishing on account of the fact that any gravitational Lagrangian is an already 
parameterised system. Notice that we have discarded a surface integral at spatial 
infinity, which reads § da k (N\jQ> k — NQ jk y). This boundary term is irrelevant for 
all class A Bianchi types. We have thus cast the action (|i~l~D into canonical form, 



S 



d x N 



M 



p ij C n hij + Q ij C n Kij - Sj c {h, K, p, Q) 



(16) 
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The associated Dirac Hamiltonian density is 

£)d = fie + m<t> kl + + e w ng ] - (17) 

Consistency of the primary constraints ([14]) when time evolution is considered yields 
the determination of the multipliers //jy and £ kl and gives rise to the unique secondary 
constraint 

X = 2p + K T kl Q lkl ^0, (18) 

which is first class — it is indeed the generator of conformal transformations. Moreover 
the spurious canonical variables Q and K can be eliminated and only the respective 
traceless parts of the corresponding tensors remain as relevant canonical variables. 

Hitherto the Hamiltonian formalism has been completely general, i.e. applicable to 
space-times without isometries; henceforth we intend to specify this setting to spatially 
homogeneous cosmologies, that is, Bianchi models, the metrics of which can be written 

ds 2 = -N 2 {t)dt 2 + h l3 (t)Ju\ (19) 

where iV is the lapse function and oj % for i = 1,2,3 are the basis dual one- forms 
(which satisfy do/ = \C\^ oji A oj k ) that characterise the isometry group of each 
Bianchi type. We adopt the usual assumptions of Hamiltonian cosmology: we choose 
the shift vector to be zero and assume that the spatial metric be diagonal and 
parameterised with Misner's variables fi(t) and P±(t). To proceed further we perform 
a canonical transformation that we have defined in a previous work |15| the virtue 



of which is to disentangle terms stemming respectively from the pure i?-squared and 
conformal variants of the general quadratic theory. The new set of canonical variables is 
{/i, P + , n^, n + , n_; Q + , <2_;P + ,'P_}, where Q± and V± parameterise the traceless 
tensors <2 T4J and Kjj respectively; the n's are conjugate to Misner's variables and the 
Vs are conjugate to the Q's. Hence the spatial integration in the action ( |16|) can be 
performed; the appropriate canonical action for Bianchi cosmologies is thus given by 



S 



J dt [n M/ i + n+/? + + + v + q + + v_ s_ - nx c - a c ^ c ] , (20) 

where the explicit forms of the first-class constraints !K C ~ (super-Hamiltonian) and 
ip c pa (conformal constraint) depend on the Bianchi type considered. For Bianchi type 
I they reduce respectively to 

"K c = --= [n + P + + n_P_ - 2Q + (Vl -V 2 )+ - ]^-^(Q 2 + + Q 2 ), (21a) 

<p c = Tl ll + T + Q + + V-Q-. {21b) 

We need to impose a gauge-fixing condition associated with the conformal constraint 
( [21 The introduction of the condition \i w as an additional constraint has the effect 
of turning the conformal constraint (|21fl ) into second class. Hence we are allowed to 



eliminate the pair of canonical variables \i and n M . On the other hand, in accordance 
with the form of the Bianchi type I metric @ we can parameterise the lapse function 
as N = e M . Therefore, we obtain the final form of the canonical action (12C 



S= dt 



U + (3 + + n_/3_ + V+Q+ + VM- - K 



(22) 
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where the super-Hamiltonian ( |21o| ) is now given by the expression 
1 



YI+V+ + Ii-V- - 2Q + (Vi -Pt)+ AV+V-Q- 



- 5 



(23) 



Varying the action ( |22|) with respect to the remaining canonical variables and their 
conjugate momenta we derive the canonical equations of conformal gravity: 



Q- 

r + 

V- 



{TI+ + AV-Q- -AV+Q, 



V+V- + Q- 



This system supplemented with the super-Hamiltonian constraint (|2 
to the system of fourth-order field equations (I 



HE 



(24a) 

(246) 
(24c) 

(24 d) 

(24e) 

(24/) 

is equivalent 
However, instead of the 



Bach equations we have now at our disposal the nice differential system (|24^) (|24?) 
and the algebraic constraint (|2~3|) . We shall seek out exact solutions by performing a 
global involution algorithm on appropriate extra constraints that yield closed constraint 
algebras. 



3.2. Global involution of extra constraints 

The involution method consists in applying the Dirac-Bergmann consistency algorithm 
to our classical system, with the Poisson bracket defined with respect to the canonical 
variables (3±, U±, Q±, V±, and after suitable conditions have been imposed. For a 
rigorous treatment on the concept of involution as applied to constrained systems, we 
refer the interested reader to Seiler and Tucker's article p3|| . Strictly speaking, the steps 
of the global involution algorithm are: 

(i) Impose an appropriate extra constraint on the canonical variables. 

(ii) Require that constraint to be preserved when time evolution is considered. This 
gives rise to secondary constraints and possibly to the determination of the Lagrange 
multiplier associated with the extra constraint. 

(iii) Repeat Step (|n|) until no new information comes out. 

Once the involution algorithm has been performed we can classify all the constraints 
into first class and second class and proceed further to the analysis of the particular 
system. 
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3.2.1. Constraint <p p ~ 0. The first extra constraint we consider expresses that the 
ratio of the variables V± be constant, namely 

I (25) 



V- - pV^ 

,(°) 



0. 



p G 



The Poisson bracket of tp p and "K c yields the secondary constraint 

2 



4 ] 



V6 



p(p 2 -3)ri-(Q--pQ, 



0. 



(26) 



The Poisson bracket of (p^ and JCd := % c + X p (f^ leads to the determination of the 
Lagrange multiplier, viz. X p « (IT_ — pU + )/[y/6(l +p 2 )], which turns out to be constant, 
and the algorithm stops. Both constraints (|25|) and (|2~6|) are second class since their 
Poisson bracket is equal to 1 +p 2 . Thus we can eliminate the corresponding spurious 
degrees of freedom. To this end we perform the canonical transformation 

Q±-* Q±±pQ T - (27) 
The action (E2T) then reduces to 



S 



dt 



where the super-Hamiltonian (|23| ) is now given by 



<K c = -—(ji + + pii_)n 



p\p 2 - 3f 
. 3(1 + p 2 ) 



V 4 



2(1 +p 2 



(28) 



-(29) 



V6 V 

Varying the action (p8|) with respect to the pair of canonical variables Q + and V + we 
obtain the canonical equations 



V- 



1 

"71 

2 /3p 2 - 1 



n++ P n_) - 
vl + 



4 /3p 2 -l', n ^ 



i + r 



4p 2 (p 2 - 3) 
. 3(1 +p 2 



2-i 



1 +' 



(30 a) 



(306) 



+ v^V l+P 2 J' + ' 1+P 2 ' 
Taking into account the explicit form of the super-Hamiltonian (p9|) we obtain a 



nonlinear first-order differential equation for P+(t), 



P 2 



(31) 



■v " ' " + v/7) v i + p 2 

which can be written as one binomial equation of Briot and Bouquet, 

u 2 = u 4 + 7 3 M , (32) 

in terms of a new variable u defined by the homography u = ±2V + ^/ (1 + p 2 )/Q, with 
7 3 = =F § (n+ +p 2 ) -1 ^ 2 . Direct calculation shows that its general solution is 

given by 

7 3 

Mt) = (33) 



4W 

where W stands for the Weierstrass elliptic function p(t — to; g^ #3), with invariants 
P2 = and #3 = — ^(Il + + pU) 2 / (1 +p 2 ), and with one arbitrary constant t (see [^j 
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pp 627 ff). In terms of the canonical pair (Q + ,V + ) the representation 
to the expressions 



Q+(t) 



V6 (n+ + P u. 

72 W 2 



1 - 3p 2 



V (t) - v / 6/n ±± pr^\ 1 



L(i +p 2 
1 

W 



dln>V 

dt 



corresponds 



(34) 



where 



dlnW 



dt 



(n+ + P n_ 
i +p 2 



To obtain the analytic form of the metric functions (3±{t) we must integrate the right- 



hand side of the expression giving V+(t). We apply the following theorem (see [£4 
p 641): If p'{zq) ^ 0, then 

dz 



2zQ{z ) + \na(z - z ) - \na(z + z ), 



p(z) - p(z ) 

with Weierstrassian functions ((z) and a(z) defined by ('(z) + p( z ) — an d v'{z) — 
a(z)((z) = respectively. To achieve our aim we take z = t and z = t z , where t z is a 
zero of the Weierstrass elliptic function, i.e. p(t z ) = 0. Indeed, we obtain 



dt 



± 



p- 



n+ + pii. 



2t({t z )+ In a{t-t z )-\na(t + t z 



(35) 



and we can write the corresponding homogeneous metrics of type I under the form 



ds = — dt + exp 



+ exp 




(36) 



+ exp 

'l+p 2 

As a particular case of the present analysis we can specialise the above solution 
(j33[) to the axisymmetric case, for which the secondary constraint ( p6| ) reduces to the 
conditions 

V- Q- 

Owing to the canonical transformation fl27|) the new variable Q_ vanishes automatically. 
Hence the super-Hamiltonian (|2^ ) becomes 

J» 1 \ Q2 



P 



E {0,±V3}. 



(37) 



2(1 +p 2 



(38) 
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and the corresponding canonical equations are 



1 /tt „ x 4 f3p 2 - ] 

-=(n++pn_)--= -f— a 



7>4 



Vl + 



Q- 



(39 a) 



(396) 



As in the general case we find out the binomial equation of Briot and Bouquet ( pi] ) the 
solution (|3~3|) of which is valid for any value of the parameter p. Hence the axisymmetric 
solution is obtained by setting p to or ±^/~3 in equation ( |3~4l ) and finally in the metric 

(ID- 



3.2.2. Constraint ip q « 0. Consider the extra constraint expressing that the ratio of 
the variables Q± be constant, namely 

Q- - <?Q+ « 0, q E R. (40) 



The Poisson bracket of (pf^ and %d := 1K C + \ q (p K q' yields the secondary constraint 



(0) 



4Q H 



(n_ - gn_ 



0. 



(41) 



V4q 2 -l)-2qV, 

To check the consistency of the involution algorithm we consider two distinct cases: 

(i) The secondary constraint (|4]) is identically satisfied if Il_ — glT + « and 
V-(q 2 — 1) — 2qV + ~ 0. Again, in order to proceed, we split up the analysis 
into two subdivisions: 

(a) If q 7^ <7, with a = ±1, the following weak equality holds: P_ ~ 2qV + /(q 2 — 1). 
The Poisson bracket of (pq and yields the constraint 

4g(l-3g 2 )(g 2 -3), 



;i-3g 2 )A g + g(g 2 -3)Q_ 



-VI 



0. 



(42) 



v/6(g 2 - 1) 

where q ^ a/v^3, for we consider nonzero canonical variables Q±. Consistency 
then leads to the determination of the Lagrange multiplier X q , viz. 
<7(<7 2 -3)l^ f 4g(g 2 -3) 



A « 



3g 2 -l 



Q- 



v/6(g 2 - 1) 



When q e {0, av^3} the Poisson bracket of ip q and J{d vanishes identically and 
the involution algorithm does not generate any new constraint. This subcase 
corresponds to a particular case of the axisymmetric solution with Il_ w qll + 
(cf. equation (]3~7j) ). For the nonaxisymmetric case, on the other hand, the next 
step in the involution algorithm yields 



n. 



4(g 2 + 1; 



16(3g 2 -l)(g 2 + r 
V6(q+lf(q-lf 



vl 



V- 



+ • 



and the last step gives 

(V-l)(ff-3) 



VQ(q 2 



V: 



At this stage no more constraints arise. Taking into account that the super- 



Hamiltonian (|3|) is weakly vanishing we find out that Il + must be equal to 
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zero. This subcase will be discussed in § |3.2.3| , as part of our discussion of the 
generic case with U± ~ 0. 
(b) If q = a, then V + is weakly vanishing. The Poisson bracket of ipq 1 ^ and IKd 
yields the constraint 

2a 

7! 

and consistency determines the Lagrange multiplier X q . The Poisson bracket 



<p® = \ q + crQ + + —V 2 « (43) 



(2) 

of (f q and 'K<£> gives the weak equality 



/4 
\ a 



-.VI )V 



The last step in the involution algorithm yields 

2^ 

At this stage no more constraints arise. The vanishing of the super-Hamiltonian 
restricts n + to be zero. As in the previous nonaxisymmetric case with 



q cr, this subcase will be discussed in § [3.2.3 



(ii) When n_ — qll + ^ 0, the Poisson bracket of tpjp and %d yields the constraint 



(P. 



(2) - V6q(q 2 -3)Q 2 + + V6(1- 3q 2 )\ q Q + + 4Q + (V + + qV-)(qV + - V-) 



q 



+n_(P+ -qV-)+ U + (qV + +V-) « 0. 

If g = ct/a/3, the involution algorithm is closed when tpf is computed, but it does 
not lead to any exact solutions since the super-Hamiltonian is not compatible with 
the constraints (fq \ with j = 1, . . . , 5. On the other hand, if q ^ o"/v^3, we have 
not been able to close the algorithm. It turns out, however, that the involution 
algorithm in this case is useless since the system under consideration does not 
produce an integration case, as it can be shown by a local study of its analytic 
structure. 



3.2.3. Constraints Ii± ~ 0. Consider now the extra constraints expressing that the 
canonical variables Il-|- be equal to zero. Contrary to the involution of the previous 
constraints (p p and <p q , the consistency algorithm here is trivial: it is unable to produce 
any new information since there are no secondary constraints. The constraints U± ~ 
are first class; we can choose their corresponding Lagrange multipliers \± to be zero. In 
that case the constraints remain weak equations to be imposed on the physical states of 
the system. Hence the appropriate system of canonical equations is the system ([24 4 )- 
(|246|), where we set to zero. Beside equation (|24a|), the relevant equations are thus 



Q + = ^=(V + Q + -V-Q-), (44a) 
Q-= -^(7>-Q+ + P+Q-), (446) 
V + = ^=(P 2 - V 2 + ) + Q + , (44c) 
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V- = -^V+T- + Q_. (44 d) 
V6 

The super-Hamiltonian ( f23| ) is now given by 

Q + (V 2 - V%) + 2V+V-Q-] - \{Q\ + Ql). (45) 



The general solution of equations (|44fl| )- (|44tj) is easy to produce under analytic form. 
We first solve equations ( |44dj and ( |44g( ) for Q + and Q_ and write down second- 
order equations for V+ and V-, namely 3V± = 4EP±, with E := V\ + V^_. Beside 
the super-Hamiltonian constraint, J{ c ~ 0, this system possesses the first integral 
V+V- — V-V+ = 5, where 8 denotes an arbitrary constant. Making use of the 
super-Hamiltonian constraint it is straightforward to produce a scalar second-order 
equation for E, namely E = 4E 2 . We integrate this last equation so as to produce 
the solution 2E = 3p(t — t Q ;0,g 3 ), with arbitrary constants to and #3 = ^-5 2 . This 
reduces the system 3V± = 4EP ± to linear differential equations of the Lame type, 
namely V± = 2p(t — t ; 0, g 3 )V±, studied exhaustively in references [^] and |26] (if we 
adopt Ince's notations, our particular case of the Lame equation is specified by h = 
and n = 1). The general solution of the system (|44a| )-( fl4"^ ) is indeed uniform. This 
confirms the single- valuedness of a possible integrability case detected in Subsection |4T3| , 
namely that the complete system with = has the Painleve property; see Subsection 
|Q| . If we introduce t z such that the transcendental equation p(t z ) = is satisfied, i.e. 
t z is a zero of the Weierstrass elliptic function, then the general solution to the system 
under study is given by the fundamental set 

P±,i = exp[-tC(t z )]^7^, V ± , 2 = exp[+tC(t z )]^-^, (46) 

a{t) a(t) 

with Weierstrassian functions ((t) and a(t) defined by ((t) + p(t) = and &(t) — 
a(t)((t) = respectively. The solutions given by the fundamental set flUp are distinct, 
provided 7^ 0, with i = 1,2,3, which is indeed the case here when g 3 7^ 0, since e i7 
with i — 1, 2, 3, are defined by 

e 1 + e 2 + e 3 = 0, 4(e 2 e 3 + e 3 e 1 + e x e 2 ) = -g 2 = 0, 4eie 2 e 3 = 5-3. 

If #3 = = 5, then ei = e 2 = e 3 = 0, and the solutions in the fundamental set are 
not distinct. This case is better understood in the light of its analytic structure. We 



shall see below that what is referred to as the singularity family /(3),± in Subsection fO 



and |Appendix B| becomes an exact two-parameter particular solution when vanish. 
With K constant, this is 

sin AK — 2 sin cos 4K ( 1 ± 2 sin 4K) 

Q+ = ~~2 (t^tof ' Q - = ~^2 Jt~ttf ' (47) 

V+ = ±- , V- = +- . 

2 t — to 2 t — to 

Denoting x '■= t — t and integrating equations ( |24a|) we obtain type I homogeneous 

metrics under the form 

ds 2 = -d X 2 + X T^K-V3cos4K dx 2 + ^ain^WSca,^^ + ^2^4X^2 (4g) 
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4. Analytic structure of the Bianchi type I system 

4-1. The Painleve strategy 

The existence of the above solutions, whether particular solutions of the general 
differential system or general solutions of specialised systems, tells nothing about the 
integrability or non-integrability of the complete system and gives no information 
whatsoever about the mere accessibility of an exact and closed-form analytic expression 
of its general solution. This is due to the fact that the global involution algorithm of 
the extra constraints, as operated above, is not related with integrability and may even 
prove to be nonexhaustive. The integrability issue will be tackled through an invariant 
investigation method of intrinsic properties of the general solution. In particular, 
the result will not depend on specific choices of the metric, within some well-defined 
equivalence class. 

The approach advocated by Painleve proceeds from the main observation, nowadays 
frequently exemplified in various areas of theoretical physics, that all analytic solutions 
encountered are single-valued or multiple-valued finite expressions (possibly intricate) 
depending on a finite number of functions — solutions of linear equations, elliptic 
functions, and the six transcendental functions systematically extracted by Painleve 
and Gambier. Therefore, the building blocks of this process are the functions, explicitly 
defined through their single-valuedness. This emphasises the relevance of the analytic 
structure of the solution, as uniformisability ensures adaptability to all possible sets 
of initial conditions. The next step of this approach deals with integrability in some 
fundamental sense. Indeed, probing the analytic structure of a system requires a global, 
as opposed to local, integrability-related property, namely the Painleve property — we 
recall that a differential system possesses this property if, and only if, its general solution 
is uniformisable or, equivalently, exhibits no movable critical singularity. At this level, 
integrable systems are defined in the sense of Painleve as those systems that possess 
the Painleve property. The Painleve property is invariant under arbitrary holomorphic 
transformations of the independent variable and arbitrary homographic transformations 
of dependent variables. In particular, our results are invariant under arbitrary analytic 
reparameterisations of time. Moreover, in keeping with the above observation, the 
requirement for global single-valuedness ought to be relaxed so as to accommodate 
the definition to integrability in the practical sense. Upon using broader classes of 
transformations, the (unavoidable) analytic structure of the solution is easily kept under 
control and either allows, or rules out, the possibility to produce a mode of representation 
of the general integral. 

Such an approach requires the analytic continuation of the solution in the complex 
domain of the independent variable and computes generic behaviours of the solution in 
a vicinity of each movable singularity. The trivial part of the Painleve method, known as 
the Painleve test, produces necessary but not sufficient conditions for a system to enjoy 
the Painleve property and requires local single-valuedness of the general solution in a 
vicinity of all possible families of movable singularities. We recall that movable essential 
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singularities are difficult to handle, since one lacks methods to write down conditions 
under which they are indeed noncritical. For a tutorial presentation of basic ideas and 
constructive algorithms aimed at the generation of integrability conditions, we refer to 
Conte [27| and Ramani, Grammaticos, and Bountis [FJ^ . 



4-2. Analytic proof of non-integrability 



The complete differential system ( [2 4 d.) - ([24 J) admits three distinct local leading 
behaviours in some vicinity of movable singularities, that is, three distinct singularity 
families, given in [Appendix B| . The first leading behaviour, denoted /(i), is valid 



in a vicinity of the movable point t = t\ and exists only when n + n_ 7^ 0. The 
existence of the second singularity family, referred to as f(2),±, valid in some vicinity 
of the movable point t = t 2 , requires Il_(2v / 3n_ ± 3II + ) 7^ 0, respectively. The third 
and last leading behaviour, denoted /(3),±j is valid in some vicinity of the movable 
point t = t 3 and is associated with another arbitrary constant parameter, denoted K 
Appendix~B| . It exists only when AK £ {0, 7r, -tt, §, — §, Tf , Tif}- The super- 



111 



Hamiltonian ( |23|) introduces no new restrictions since, at leading order in % := t — ij, 
for all i = 1,2,3, it always vanishes in some vicinity of movable singularities around 
which leading behaviours /(i), f(2),±, and /(3),± hold. Around each species of movable 
singularities one must inquire whether it is possible to generate single-valued generic 
local representations of the general solution. 

Such local expansions will be required to be both generic and single- valued. As it is 
well known, both requirements are associated with Fuchsian indices of each family. The 
indices have been computed for each leading behaviour and all preliminary conditions 
have proved to be fulfilled (for a thorough account an interested reader may consult 
[Appendix B|) . However, near both movable points t 2 and t 3 in C, it is not possible 
to build single-valued generic local expansions, unless some constraints are imposed. 
Introducing movable logarithmic terms in these local series one may regain genericness — 
at the expense of losing the Painleve property. The complete proof of these claims can 
be found in |Appendix B . 



Such multiple-valuedness, betrayed by local investigations, also pertains to the 
general solution itself. This analytically proves that the system under consideration is 
not integrable. Its general solution exhibits, in complex time, an infinite number of 
logarithmic transcendental essential movable singularities; in other words, an analytic 
structure not compatible with integrability in the practical sense: the quest for generic, 
exact and closed-form analytic expressions of the solution is hopeless. This result holds 
under space-time transformations within the equivalence class of the Painleve property; 
see Conte p7| . Yet, local information produced by the Painleve test may still be used 
in order to extract all particular systems which may prove integrable. In the case under 
study, see Subsection |3.2| , it turns out that all particular solutions are meromorphic, but 
this is only a result of the application of the method. 
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4-3. All integrable particular cases 

Even though this does not provide a proof but merely an indication of uniformisability, 
the Painleve test requires local single- valuedness around all possible species of movable 
singularities. We shall explicitly use the results given in |Appendix B| . Near the movable 



point t = t 2 , the complete set of integrability conditions is satisfied if one imposes 

n_ = ±Vm+ (49) 

and requires the vanishing of the arbitrary constant parameter associated with index 
j = —3, i.e. Q+ -3 = 0. Both conditions must be required simultaneously. On the other 
hand, if one sets 2v / 3fl_ = =F3I1 + or EL = 0, the family f(2),± dies out and ipso facto 
no restriction needs fulfilment. 

Near the movable point t — t 3 , single- valuedness may only be recovered in two 
particular cases. The first case deals with a local representation of the general solution of 
that specialised system obtained with U± = 0. In this case, local leading behaviours /(i) 
and f(2),± die out whilst, near t = ts, a meromorphic local generic expansion is produced. 



Indeed, in this case, § |3.2.3| has already provided the uniform general solution in closed 
analytic form. The second case deals with a local representation of some particular 
solution of the complete differential system and requires that the arbitrary parameter 
K be such that 

sin4if=^ n+ cosAK = ± -. U ~ (50) 

5. Conformal relationship with Einstein spaces 

A four-dimensional space-time (M, g) can be mapped onto an Einstein space, under a 
conformal transformation, g ab = e 2a ^g a b, if and only if there exists a smooth function 
cr(x) that satisfies 



1 _ A 

J ab 



V a( rV t a - V Q V 6 <7 - -ga b g mn ^mO-V n (r - -e 2a g ab , (51) 

z o 

where the tensor L ab is defined by equation (|5|) and A := R/4 denotes the cosmological 
constant characterising the Einstein space (M, g) fl29|. The first integrability condition 



of equation ( pll ) is given by 

V m C m abc + C m abc V m a = 0. (52) 

This is the necessary and sufficient condition for a space to be conformally related 
to a "C-space" ||3Cf| . The second integrability condition of equation (|5l|) is nothing 
else than equation (§): fulfilment of the Bach equations is a necessary condition for 
a space to be conformally related to an Einstein space. If considered separately the 



above integrability conditions of equation ( pl| ) are merely necessary conditions with 
regard to the conformal relationship with Einstein spaces. However, Kozameh et al. 
have proven that they constitute a set of sufficient conditions as well: A space-time 
(M, g) is conformally related to an Einstein space (M, g) if and only if equations (0) and 
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are fulfilled |3T| . 



Now let us specify the components of equation ([51]) for the Bianchi type I metric 
(|); we obtain the differential system 

A 



To - a 2 + 5(/£ + Pi) - 



0. 



a 2 + 2&(p+ + V3$_) +/3 + + -$l-$i 



A 

3 
A 



e 2CT = 0, 



a 2 + 2&0 + - y/3$_) -p + - -ft -Pi- -e 2CT = 0, 



A&/3+ - 2(3+ -Pi -Pi- 



2rr 



0. 



(53 a) 
(536) 
(53c) 
(53 d) 



which readily integrates to yield the necessary and sufficient condition for a space of 
Bianchi type I to be conformally related to an Einstein space, namely 



V± = -V6p± = -V6k±e 



2a 



(54) 



and provides us with a differential equation that enables to determine explicitly the 
conformal factor x := e 2cr , viz. 



x 



-Ax 3 + A(k 2 + + k 



(55) 



This last equation can be written as the Weierstrass ODE in terms of the Weierstrass 
elliptic function W := p(t — t ; g 2 , #3), with invariants g 2 = and #3 = — |A 2 (fc^_ + kl). 
Therefore, the conformal factor that brings a Bianchi type I space onto an Einstein 
space is given by 

Taking this last result into account we get from equation (|54|) the explicit form of the 
functions V±(t), that is 



P±(t) = -^Ak ± W-\ 



(57) 



So far we have not used the Bach equations ( [A.1|) ~ ([A.4|) nor the equivalent canonical 
system ( p4c| ) (|24^ ) . If we do so, we see that the expression (|57|) does actually coincide 
with the solution ([34]) upon identifying 



_ A; + n + + k_U_ 

V= k~+' A{k 2 + + k 2 _) ' 

Thus we can rewrite type I homogeneous metrics (^) under the equivalent form 



(58) 



ds 2 = -dt 2 



+ exp 
+ exp 
+ exp 



2{k+ + V3k^ 
^/kl + k 2 

2(k+ - y/Zk- 



Ak+ 



k 2 

Kit 




(59) 
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This is in agreement with the fact that imposing the constraint (|25|) — here, a direct 
consequence of equation ( |54j ) — on the Bach equations is equivalent to requiring condition 
(|52|) to be fulfilled, as it can be proved with the help of Reduce. In other words, 
assuming that the ratio of variables V± is constant is a necessary condition for a Bianchi 
type I space to be conformally related to an Einstein space; in conformal gravity, it 
becomes a sufficient condition as well. Thus the only way to find out a solution to 
the Bach equations that is not conformally related to an Einstein space is to relax the 
constraint (|25|) . In accordance with our analysis of the preceding sections, we have 
indeed obtained the only closed-form analytical solution that cannot be mapped onto 
an Einstein space, namely the general solution fl46[) to the canonical system ( j!4dV(|444) , 
thereby confirming Schmidt's conjecture on the existence of such solutions ||32| . 

Now let us examine the particular case of a vanishing cosmological constant A. 
Equation (|55|) becomes trivial and yields (up to an irrelevant constant of integration) 

e 2 ^=t-t =: X . (60) 

Inserting this result into equation ([54]) and integrating we obtain type I homogeneous 
metrics under the form 

ds 2 = ~d X 2 + X 2(fc++v^-) dx 2 + X 2(k + -V3k.) dy 2 + x ~^dz 2 , (61) 



which coincides with the exact 2-parameter particular solution fl48|) in the case of zero 
n ± , upon identifying k + = =f^ sin 4_K" and k_ = — |cos4i^. We then perform a 
conformal transformation with conformal factor (|60|) and introduce the proper time 
T ._ t^ 3 / 2 ; hence the metric (|61] ) that characterises the conformal Einstein space with 
zero cosmological constant is simply 

ds 2 = -dT 2 + r §[2(*+W3M+l] dx 2 + r §[2( fe+ -v^-) + l] d2/ 2 + T l(l-4k + ) dz 2_ ^ 

If we denote by p i7 with i = 1,2,3, the successive exponents of r in the spatial line 
element of the metric (|6~2}), we obtain 

1 



Pl+P2+P3=l, PI+PI+PI ,;, 



i + KK + ki) 



(63) 



Computing the Bach equations for the metric (|6l|) we derive a constraint on the 
parameters k±, namely k\ + k 2 _ = j (and thus p 2 + pi + p 2 = 1), that leads to the 
vacuum Bianchi type I Kasner solution in general relativity. 
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Appendix A. Bach tensor of the Bianchi type I space-time 

We have computed the components of the Bach tensor (Q) for the Bianchi type I metric 
(0), with the ExCALC package in Reduce. Up to irrelevant multiplicative factors, they 
are: 

B 00 = 2'(3_P_-(3 2 _- 12/31 - 2AP 2 J 2 + + 2 p + 12#, (A.l) 
B n = \/3(J3_ - 24/L/?! - 8/L/3+ - 16/L/3+/3+) - 2 "p_ /3_ + (3 2 _ 
— 16/3_/3_/3 + + 12/31 - 8/3^/3+ + 24/3^ + 0+ 

-2 fi + P + + P 2 + - 24/3 + /3* + 12/3 + , (A.2) 
B 22 = VS0- - 24/L/?! - 8/L/3J - 16/L/3+/3+) + 2 /?_/?_- /3! 

+16/3_/3_/3+ - 12/31 + 8/3^/3+ - 24$ 2 _$ 2 + - 0+ 

+2 p + P + -(3 2 + + 24P+P* - 12/3 + , (A.3) 
B 33 = 2 /3_ 0_ - /3! + 32/3_/3_/3+ - 12/31 - 16/3^/3+ - 24/3^/3^. 

+2/3+ + 2 /3 + $ + - ft - 48/3 + /3 + - 12/3 + . (A.4) 

Appendix B. Painleve analysis 

Throughout this appendix we shall consider the auxiliary variables 

N {n-N)j-+j+ 

n=0 j=nj— 
N (n-N)j-+j+ 

v ± =v£= y £f Yl p{ ±l * j+P± ■ ( B - 2 ) 

n=0 j=nj— 

Local analytic results presented here have been obtained through Reduce [|18| 
implementations of the Painleve test ||3 3) | . The complete differential system possesses 
three distinct singularity families, referred to as /m, f(2),±, and /(3),±- 111 the series (P-l| ) 



and (|B.2|) we have introduced: 

• the integer exponents a± and b±, taking their values amongst { — 1,1} and inducing 
homographic transformations of dependent variables Q± and V± (indeed, the 
Painleve property is known to be preserved under such mappings); 

• the order N of the perturbative Painleve test, taking its value amongst {0, 00}; 

• whenever N 7^ 0, the small perturbative complex parameter e; 

• the generically complex numbers j_ and j + , denoting respectively the lowest 
and highest Fuchsian indices — Fuchsian indices give the orders at which arbitrary 
constant coefficients enter the double expansion ( |B.1| )-( |R2|) ; 

• Q±j and V±j, the constant coefficients of the double expansion; 

• the expansion variable x °f the inner Laurent series, given by t — t\, t — t 2 , and 

t-t 3 m the cases f w , /( 2 ),±, and /( 3 ) )± respectively; 
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Table Bl. Families of movable singularities for the Bianchi type I space-time in 
conformal gravity. 



SF 


a±,b± 


N q±,p± 


-6(°) 

«±,0' ' ±,0 


Indices 




1(1) 


a± = — 1 
6± = -l 


oo q± = — 1 
P± = -2 


n(°) - \/6 

-5(0) 2^6 


{-2,-2 


;-i s -i} 

= (-2,-1) 



/(2),± a+ = -1 


oo q+ = 


-l 


qV 0) o 


i_3%/2 

— 31 n_ 


{-3;-l;l;3} 


a_ = -1 


q- = 


-l 




_ 9\/2 

2y^n_±3n + 


(.7-,i+) = (-3,3) 


6+ = l 


p+ = 


-l 




_ 

4 




= 1 


p- = 


-l 




_ x 3v^ 
r 4 





/ (3) ,± a+ = l g+ = -2 efj, = -#(l±s-2 S 2 ) {-1;0;3;4} 

a_ = l q_=-2 Q (0 . ) = -^c(l±2 S ) (j'- J J+) = (-l J 4) 



6 + = l p+ = -l PS=±^ 

6_ = 1 p_ = -1 ^i 0) = 



the singularity orders, i.e. leading powers g± and p±; 
the nonvanishing coefficients Q^ and V± - 



For each family, all aforementioned parameters are given by entries in table pi] , 
where SF stands for 'singularity family' and with s := sin4i^ and c := cos4.fr. In the 
case /(3),± the reason for the peculiar parameterisation of leading coefficients is that it 
produces transparent integrability conditions; see below. The rightmost column provides 
the set of Fuchsian indices. 

Preliminary integrability conditions are fulfilled near all movable singularities, 
located at t = ti, for all i = 1,2,3. Indeed, singularity orders, i.e. leading powers q± 
and p±, are negative integers; Fuchsian indices are integer numbers; maximal families 
are generated, thus producing local representations of the general solution, and rank 
conditions are satisfied by each set of Fuchsian indices — the multiplicity of each distinct 
Fuchsian index equals the number of arbitrary constant parameters to be introduced in 
the expansion at current Fuchsian index. 

Perturbative approaches to the Painleve test are required: the approach due to 
Ablowitz, Ramani, and Segur ||34| , ^], in the case f(3),±, and the Fuchsian perturbative 



approach devised by Conte, Fordy, and Pickering ||27], |35|, |3q|, in the cases /(i) and /(2),±- 
The ansatz previously introduced enables us to generate some integrability conditions: 
for each value of (n, j) such that ^ n ^ iV and j is an index, an orthogonality condition 
must be fulfilled. Our results may be summarised as follows. 

For all ^ n ^ 5 and all j G {—2,-1}, the integrability conditions are fulfilled 
near t — t\. The family fci) does seem to generate a single- valued local representation 
of the general solution, but we recall that the Fuchsian perturbative method is not a 
bounded algorithm — stricto sensu this is not an algorithm. However that may be, this 
leading behaviour does not seem to be informative from the integrability point of view. 
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On the contrary, dealing with the family /(2),±, a local expansion near t = t 2 requires 
the presence of logarithmic terms in order to re-establish genericness. The (n, j) = (0, 3) 
integrability condition is violated, unless one imposes the constraint 

= 0. (B.3) 

Moreover, the constraint Il_ = ±y/3U + does not preclude a subsequent incompatibility, 
encountered when (n, j) = (1,3). This last incompatibility vanishes only if one freezes 
the arbitrary constant parameter associated with the index j = —3 so as to satisfy the 
relation 

^S«snl = * Q« 3 = 0. (B.4) 

Since the multiplicity of index j = — 1 is one, it is always possible upon a gauge choice 

to impose one scalar constraint amongst the two arbitrary parameters associated with 

that index j = —1, namely, t 2 and Q+Li- We shall use this gauge freedom and set 

Q+ _i = 0. Under these circumstances, the Fuchsian perturbative method reduces to 

the test of Ablowitz, Ramani, and Segur — the algorithm stops. 

We now consider the family /(3),± and produce integrability conditions in a vicinity 

of the movable point t — t%. Again, logarithmic terms are required if one demands 

generic representations. The (n,j) = (0,3) integrability condition is incompatible, 

unless the constraint 

n + cos AK =F n_ sin AK n . , 

= (B.5 

sin AK cos AK 

is obeyed. The last occurrence of an index, namely (n,j) = (0,4), gives no new 

information — the index j = A turns out to be compatible. Local single-valuedness 
may only be recovered in two particular cases: 

(i) In the cases where both n + and Il_ vanish the general 4-parameter solution is 
locally single- valued near the movable point t = i 3 . Indeed, when ILj- = the only 
surviving family is /(3),± an d the weighted-homogeneous leading behaviour becomes 
an exact two-parameter meromorphic particular solution; see § 3.2.3 . 

(ii) When Il + ^ and Il_ ^ 0, local single-valuedness may only be recovered in the 
particular case 

sin4if := — n+ cosAK:=± - U ~ (B.6) 

Vn 2 + + iP_' V n + + n - 



Integrability results are summarised in table which gives, for each singularity 
family: all relevant existence conditions; nontrivial integrability conditions — i.e. incom- 
patibilities of conditions for the system to possess the Painleve property — ; and all 
possibly integrable particular cases. The column providing the set of incompatibilities 
gives the value of (n, j) where nontrivial conditions were met. 

Near the movable point t = t 2 , the super-Hamiltonian constraint CK C ~ requires 
2+i = 0- Taking into account all /(2),± integrability conditions one thus generates 
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Table B2. Local integrability conditions for the Bianchi type I space-time in conformal 
gravity. 



or 


J_J-A.lSliCllL.C L.U 11 LlllilU IIS 




dons 


Integr 


Q K O PQOOO 


/(i) 


n+n_ ^ o 











/(2),± 


n_(2V3n_ ±3n + ) ^ o 


and 


= (0,3) 
= (1,3) 


n_ = 

and 


±V3n+ 
= o 


/(3),± 


4K^{0,7r,-7r,f,-f,Tf,Tf } 




= (0,3) 
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two-parameter meromorphic representations that begin like the truncated sums 
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Remarkably enough, if we set Q?^ 



^ in these expansions, we are brought back to 
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manifolds with axial symmetry, since truncated-from-below Laurent series for V- and 
Q + then identically reduce to the expansions for ±y/3V + and iy^Q- respectively. 

Near the movable point t — t^, the super-Hamiltonian constraint "K c ~ requires 
Pi°4 = 0. Let us denote cr := + H 2 _ and consider the family /(3),± in the integrable 
case such that a/it sin4.fr = il + and t/ctcosAK = ±n_. In such cases, the local Painleve 
test produces the two-parameter meromorphic series 
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Again, we note that if one imposes 8V_ 3 :— ±\/2II + in these expansions, local 
representations of axially symmetric manifolds are produced. 

To end this appendix, let us note that in cases such that both 11+ and Il_ vanish, a 
fourth singularity family, hereafter referred to as fu),±, may enrich the analytic structure 
of the system. In terms of the expansion variable x '■= t — t^, this fourth local behaviour 
has leading terms 

S + ~±^S<V> 8-»SV, V,~~X-\ P-»T^V\ (B.7) 

and four Fuchsian indices, namely {— 4;— 1;0;3}, for which the rank integrability 
condition is satisfied. However, when the super-Hamiltonian constraint "K c ~ is 
to be satisfied, this implies the restriction =f6\/2Q_o = ^ an< ^ the family /(4),± dies 
out. Therefore, the new family enriches the analytic structure of the system only when 
% c ^ 0. 
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